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LASER- VELOCIMETER SURVEYS OF MERGING VORTICES 


IN A WIND TUNNEL - COMPLETE DATA AND ANALYSIS 
Victor R. Corsiglia, James D. Iversen,* and Kenneth L. Orloff 

Ames Research Center 


SUMMARY 


The merger of two corotating vortices was studied with a laser velocim- 
eter designed to measure the two cross-stream components of velocity. Mea- 
surements were made at several downstream distances in the vortex wake shed 
by two semispan wings mounted on the wind-tunnel walls. The velocity data 
provided well-defined contours of crossflow velocity, stream function, and 
vorticity for a variety of test conditions. Downstream of the merger point, 
the vorticity was found to be independent of the downstream distance for 
radii smaller than r/b = 0.05. For larger radii, the vorticity depended on 
the distance from the wing. Upstream of the merger point, a multicell vorti- 
city pattern was found. 


INTRODUCTION 


In recent years, considerable research into means of deintensifying the 
wake of large aircraft has taken place (ref. 1). Some of this work has 
focused on the process whereby corotational vortices merge as a technique for 
diffusing the wake vorticity. Experiments in this area have shown (refs. 2-4) 
that when the trailing-edge flaps of a large transport aircraft were set so as 
to shed multiple vortices per side, the vortices can merge to form a diffuse 
wake so that reduced rolling moments are imposed on a following aircraft. 

Rossow further has found (ref. 5) that a fin mounted on a wing upper surface 
so as to introduce an additional vortex per side into the wake also results 
in a greatly diffused wake. In addition, his scheme allows the wing trailing- 
edge flaps to be set at their usual landing position. 

Several methods have been used to investigate the merger process theoreti- 
cally. The point vortex Interaction approach provides simplicity and contains 
the essential ingredients of convection. Rossow (refs. 6, 7) and Christiansen 
(refs. 8-10) have used this inviscid approach to show that a boundary exists 
between merger and nonmerger, depending on the relative strength, size, and 
separation of the two vortices. Another method, second-order closure theory 
for turbulence, provides a greater sophistication in the modeling at the 
expense of Increased computational complexity. Using such a code, Bllanin 
(refs. 11 and 12) has shown that maximum turbulent kinetic energy is higher an 
and maximum vorticity is lower during the merging of vortices than during the 
decay of a single isolated vortex. He noted that the presence of two vortices 
destroys the axial symmetry and allows turbulence to be produced. 

*Professor, Iowa State University, Ames, lA. 



An intermediate technique is that of Steger and Kutler (ref. 13) wherein 
diffusion is included by use of a simplified viscosity model. 

A need exists for detailed data describing the merger process in order to 
assess the relative merits of the various computational schemes. In one set 
of experiments directed at this need, merger was studied by flow visualization 
(refs. 14-16) using the same wind tunnel and models as were used in the present 
study. Two semispan wings mounted on the wind-tunnel walls provided two 
vortices whose interaction and merger were observed. The distance from the 
wings to the merger point correlated well with the initial vortex spacing, 
suitably normalized by a vortex shape parameter termed the circulation defect 
diameter. No boundary between merger and nonmerger was found, suggesting a 
significant influence of turbulent diffusion. However, the experimental data 
did not greatly exceed the critical spacing predicted by Rossow (ref. 6) and 
Christiansen (refs. 8-10). 


The objective of the present study is to provide additional experimental 
data for various vortex merger processes in order to compare and evaluate the 
theoretical models described above. Time-averaged velocities were measured 
using a laser velocimeter on crossflow planes at several distances downstream 
from the wings. As a first step in the evaluation of the theories, a companion 
paper (ref. 17) uses the intermediate technique of Steger and Kutler (ref. 13) 
to calculate results for comparison with the present data. Additional results 
of the present study are contained in references 18 and 19. 

The assistance of Professor P. Raj and student R. Brickman in recording 
the measurements and reducing the data is gratefully acknowledged. 


WIND-TUNNEL FACILITY AND LASER VELOCIMETER 


The data presented herein were obtained with a laser velocimeter 
installed at the Iowa State University 76- by 91-cm (30- by 36- in.) Open 
Circuit Wind Tunnel. This tunnel incorporates an extended test section 
length of 6 m in which the turbulence intensity was measured 0.1% over most 
of the test section with some local regions as high as 0.9%. The flow fields 
of interest were generated by two semispan wings mounted on the bottom and 
side walls of the test section at any one of four streamwise locations 
(fig. 1). These wings were identical NACA 00] 2 airfoils with constant chord 
and elliptical tips. The distance that each wing protruded into the test sec- 
tion was adjustable so that the relative position of the shed vorticity could 
be varied. The Installation of the wings was identical to that of earlier 
studies (refs. 14, 15) and additional details are contained there. 

The laser velocimeter, shown schematically and photographically in 
figure 2, was especially designed for the present experiment. It incorporated 
two dual-scatter crossed-beam systems that operated independently at wave- 
lengths of 488.0 (blue) and 514.5 (green) nanometers originating from a single 
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argon-ion laser. This instrument differed from an earlier version of the two- 
color laser velocimeter (ref. 19) only in the arrangement of the optical com- 
ponents, the use of fixed focal length optics rather than scanning optics, and 
the use of acousto-optic cells in each channel to eliminate directional ambi- 
guity. Inclining the green output channel at 30“ to the blue in the vertical 
plane arranged the optical components so as to make the velocimeter sensitive 
to the two crossflow components of velocity. 

The lateral v and vertical w components of velocity were obtained 
by the following combination of the velocities measured by the blue and green 
channels of the velocimeter. 

V = -2 [vq cos 6 - Vg cos (30 - 0)] 
w = 2[vq sin 9 + Vg sin(30 - 6)] 

where 6 in the inclination of the centerline of the blue optics channel 
(9 = 10“). The streamwise velocity u was measured by rotating either beam- 
splitter by 90“ to align the plane of the output beams in the x-direction. 

The amplified outputs of the photomultiplier (PM) tubes (fig. 3) were 
analyzed with a two-channel frequency tracker described elsewhere (ref. 20). 
The tracker was used to demodulate the velocimeter signals, thus providing 
voltages proportional to the PM tube signal frequencies. The voltages were 
then passed to an analog scaling computer that was used to resolve the above 
equations into the w and v velocity components. 

Using the method described in reference 21, the spatial resolution of the 
optical system was estimated by computing the probe volume size and accounting 
for the off-axis receiving optics. The dimensions of this sensing volume were 
computed to be 0.3 x 0.3 x 1.8 mm. Accuracy of the mean velocity measured 
within this probe volume was determined from the frequency tracker and analog 
scaler specifications to be within 2% for the v and w components of 
velocity. 

Vertical translation of the probe volume was accomplished by raising and 
lowering the entire velocimeter on a hydraulic table. Lateral motion was 
accomplished by manual positioning of the system on rails with ball-bushing 
support. Positions were reproducible to within 1.0 mm. 

To ensure the acquisition of nearly continuous signals, an aerosol of 
mineral oil was injected upstream of the wind-tunnel inlet. By strategically 
placing the aerosol outlet the region of interest could be adequately seeded. 


PROCEDURES AND CONDITIONS 


For each test, both wings were set at the same span and angle of attack 
(a = 10“), and the wind-tunnel speed was set to approximately 46 m/sec 
(Reynolds number = 4.7x10^ based on chord). 
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In earlier studies (refs. 22, 23) where the velocities in the vortices 
were also measured, a high-speed traverse was required to circumvent the 
problem of meander of the vortex center due to large-scale wind-tunnel turbu- 
lence. In the present tests, the vortex meander was greatly reduced so there 
was no need to rapidly scan the region of interest. Data were taken only when 
the vortex center was observed to be in a prescribed location. To accomplish 
this, the vortex was observed with the aid of a light slit and two crossing 
lines made up of a shadow and a laser beam (fig. 4). Data were continuously 
averaged to remove high-frequency noise (time constant =0.3 sec) and sampled 
several times at each point. In this way, repeatability of the results was 
generally within 4%. Measurements were taken in 10-mm increments vertically 
and laterally at those locations where there was sufficient aerosol to seed 
the flow (20-mm grid increments for the streamwise velocity component). Up 
to 380 locations were measured at a streamwise location for a single 
configuration. 


RESULTS 


Figure 5 shows the locations at which data were obtained for the case of 
the greatest distance downstream of the merger point (x^ = 3.0). Merger dis- 
tance was calculated using the expression found empirically in references 14 
and 15. The lateral distribution of the vertical and lateral component of 
velocity is also shown for two values of vertical position. The scan through 
the vortex center shows points missing in the center due to lack of seeding 
material there. Two rows away from the center, however, there was adequate 
seeding of the flow for the entire scan. Tunnel-wall corrections were not 
applied to the data because an estimate of these corrections indicated that 
they were small compared to the measured cross-stream velocities. 

The streamwise component of velocity was measured in a separate survey 
from the crossflow velocity measurements because the laser veloclmeter could 
determine only two components of velocity at a time. The streamwise velocity 
component is presented in figures 6 and 7 for several configurations. A 
variation of up to 5% was found (fig. 7). However, this variation was not 
attributed to the vortices, since a similar variation was found when the wings 
were set to zero angle of attack (fig. 6). This result is in agreement with 
earlier studies (refs. 22-24) wherein the streamwise component of velocity 
differed appreciably from free stream only over a limited radial extent from 
the vortex axis, and that variation was associated with loss of total head 
due to friction on the wing. Ih the present study, small values of radius 
could not be attained, and negligible total head loss was expected since the 
wings were thin and smooth. 

The crossflow velocity magnitude was computed according to the expression. 

V = (v^ + w^) 

These results were plotted in bo_th the vertical and lateral directions, and 
contours of constant values of V were obtained (fig. 8). The symbols on 
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the figure indicate the points available from cross-plot. As can be seen, 
the contours are well defined over much of the region and vortex merger 
seems complete. 

The stream function shown on the lower part of the figure was computed 
as follows. For two-dimensional flow in the y-z plane, the stream function 
is defined as 

'I' “ = I -V dz + w dy (1) 

2 1 J 

1 

The integration was accomplished by fitting four consecutive velocity points 
in the y or z direction with a third-order Lagrange polynomial and inter- 
polating for a Gauss-Legendre quadrature. The integral would be exact if the 
velocity distribution were a polynomial of degree 3 or less. Also, if the 
flow were not two-dimensional, equation (1) would lead to some error since, 
from Stokes theorem and continuity, the closed contour line integral of 
equation (1) is 




-V dz + w 


dy 



dy dz 


( 2 ) 


If the streamwise velocity is independent of x, equation (2) is identically 
zero. This calculation was performed for each 10-mm square for which velocity 
data were obtained. The values for the contour line integral calculated in 
this way were not only small but also random in sign. The flow was thus 
assumed to be sufficiently two-dimensional that the stream-surface cross 
sections calculated by equation (1) were independent of the path used. 

To calculate the stream function at any one point in the flow field, a 
reference point was selected at which the value of the stream function was 
arbitrarily set equal to zero. The stream function at every point in the flow 
field was then obtained by calculating the change in stream function according 
to equation (1) along the combination of lateral and vertical distance incre- 
ments Av and Az which provided the shortest integration length. This has 
the effect of minimizing the error due to axial velocity gradient (eq, (2)), 

As can be seen, very smooth, nearly circular, stream functJ.on contours were 
obtained. The smoothness of the plot as compared to the V plot reflects the 
smoothing effect of the calculation. 

Figure 8(c) shows contours of constant vorticity. To obtain these values 
the data were differentiated using the following four-point expression: 

3w _ ^ 

^ 3y 3z 

= -^ [ (W2 - wi + W3 - W4) - (V4 - vi + V3 - V2)] 
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where the subscripts refer to the four corners of a measured grid element. 

These results show (fig. 8(c)) a lower degree of axial symmetry in the vorti- 
city than in the stream function or velocity. An interesting feature is the 
wave pattern (approximately seven periods per revolution) around the vortex. 

A similar feature was noticed by Mertaugh et al. (ref. 25) in their hot-wire 
anemometer surveys of the near-field wa'Ae of a full-scale aircraft in flight. 

A calculation was made to determine the magnitude of the error that might 
occur in the numerical calculation of vorticity. For this purpose, an equation 
for the circumferential velocity was chosen that approximates the velocity 
values shown in figure 8: 


= -1.22r + 0.25 (3) 

From this expression, both vertical and lateral components of velocity were 
computed at a set of grid points to simulate the measured data. The vorticity 
was then computed numerically by use of the procedure given earlier in the 
paper and by use of an analytical expression derived from equation (3): 

^c 

r dr 

These two vorticity results are compared in figure 9. As can be seen, the 
numerically computed values differ only slightly from the exact circular con- 
tours. It was, thus, concluded that the numerical procedure was accurate 
enough for the present purposes. 

Figure 10 shows a set of contours of crossflow velocity, stream function, 
and vorticity at a streamwise location nearer to but still downstream of the 
merger point ~ 1*2. Although the stream function contours appear smooth 
and circular, the velocity and vorticity both lack axial symmetry when compared 
with the previous figure. In fact, a significant second cell of vorticity is 
present at = 1.2. Merger has been defined as the coalescence of vortex 
centers based on flow visualization. If merger were defined as the coalescence 
of centers of vorticity, then the vortices shown in figure 10 would not be 
merged. 

The comparable results at a streamwise location even closer to the merger 
point (x^ = 0.8) are presented in figure 11. To obtain data at this value of 
X , the span of the wings was reduced by 3% to achieve an increase in the 
initial vortex spacing, and the distance from the wings to the measurement 
point was increased to 7.4 spans. The shapes of the velocity and stream 
function contours are similar to those in figure 10. The wave pattern around 
the vortex in the vorticity contour is still evident. 

To compare the results shown in figures 8, 10, and 11 with each other, an 
average radius for each contour level of the vorticity was obtained. This was 
done by equating the area within each contour level (ignoring the secondary 
cells) to the area of a circle. The distribution of vorticity with average 
radius forms a single curve for all three configurations (fig. 12). These 


6 


results show that within f = 0.05 there is no evidence of diffusion of 
vorticity over the range of distances tested. A Lamb vortex profile fits this 
curve well to a radius Y = 0.02: this profile does not strictly apply, how- 

ever, because it represents the decay of a line vortex instead of the decay 
of two distributed vortices in proximity of one another, which is the present 
case . 


A similar comparison can be made by computing the circulation. To obtain 
these values, the component of velocity which is tangent to a circle of radius 
r" was computed according to the expression 


V 


c 


/wy - vz V^^ 
\y2 + z’2/ 


The area within contours of circumferential velocity (not shown) was then 
measured to obtain the average radius from which the circulation was computed 
as follows: 


r = 2TrrVj, 

The results show (fig. 13(a)) that, to a radius r = 0.05, the circulation 
profiles are also nearly identical. This agrees with the result shown on the 
(frevious figure. However, at larger radii there is evidence of outward diffu- 
sion or convection of vorticity as seen by the lower circulation values. The 
results for 3c =7.2 and 7.4 nearly superimpose and are lower than the results 
for 3c =5.8. On the other hand, the distances from the merger point are not 
in the same order as the circulation curves. It can be concluded that the vor- 
tex decay process depends primarily on the distance from the wings. 

For comparison, the circulation profiles for a single wing are shown 
(fig. 13(b)). The circulation rises slightly with downstream distance indi- 
cating that the rollup process is still occurring. Also shown is a calculated 
profile using the method of Betz as discussed by Rossow (ref. 26). The Betz 
profile agrees with the measured profile for the maximum downstream distance 
5^ = 7.2. Since no decrease in circulation with increasing downstream dis- 
tance was found for the single-wing case (in contrast with the decrease found 
for the two-wing case), the vortex decay shown in figure 13(a) must be the 
result of the presence of the second vortex. 

Figure 14 shows the results at a streamwise location upstream of the 
merger point = -4.2. The separation of the two vortices is clear and 
agrees with the location obtained by flow visualization. The crossflow veloc- 
ity vanishes in the region midway between the vortex centers as would be 
expected. An Interesting feature of the data presented in figure 14(c) is the 
multicell pattern of the vorticity contours, possibly caused by incomplete 
rollup of the vorticity from each wing, since the single wing results showed 
incomplete rollup at = 3.0 (fig. 13) compared to 3c^^ = 1.8 in figure 14. 
Since the distance between the main vortices in the pattern is comparable to 
the distance between the small vortices and their nearest neighbor, the merger 
process and the initial rollup process appear to be occurring simultaneously. 
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ADDITIONAL CONFIGURATIONS 


Contours of velocity, magnitude, and vorticity for those configurations 
in table 1 not presented earlier are contained in figures 15-17. 


CONCLUSIONS 


Detailed surveys were performed of the velocity field of two vortices 
undergoing a merger process. Sufficient velocity data were obtained to define 
the crossflow stream function and vorticity over a plane at several streamwise 
locations both upstream and downstream from the merger point. The axial sym- 
metry of the merged vortex decreased as the merger point was approached from 
downstream with secondary cells evident at the nearer locations. Downstream 
of the merger point within a small radius from the vortex center, the vorti- 
city was shown to be independent of streamwise location. At larger radii, 
the reduction in vorticity with downstream distance was found to depend on 
the distance from the wing rather than the distance from the merger point. 

A single isolated vortex showed no aging over the streamwise distances 
studied. It was therefore concluded that the second vortex was responsible ^ 
for the reduction in vorticity. At a streamwise location forward of the 
merger point, a multicell vorticity pattern was found. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, June 21, 1978 
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TABLE 1.- LIST OF CONFIGURATIONS 
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(b) Photograph. 
Figure 1.- Concluded. 
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(a) Schematic of optical arrangement 
Figure 2.- Laser velocimeter. 



(b) Photograph of laser veloclmeter installation. 
Figure 2.- Concluded. 
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Figure 3.- Schematic of the electronic signal processing 













Figure 4.- Light slit photograph to illustrate flow visualization of the 

vortex center. The vortex center appears dark due to the absence of seed 
ing material. 
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Figure 5.- Grid of data obtained and typical measured results of velocity 
survey; b = 749 mm, = 7.2, = 3.0. 






y, mm 


Figure 6.- Streamwise velocity component with no vortex in the ’ 
two horizontal traverses; y and z reference same as in figu 
a = 0°, x^ = 3.2, b = 72.4 cm, 2 wings. 
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(a) ^ = 7.4 

Figure 7.- Streamwise velocity component with vortex in tunnel for various 
horizontal traverses at different distances from the vortex center; 
a = 10°, b = 72.4, 2 wings. 







(c) = 1.8 

Figure 7.- Concluded. 
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(a) Crossflow velocity. 



(b) Stream function. 

Figure 8.- Results calculated from the measured velocities; b = 749 mm 
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Figure 9.- Comparison of vorticity magnitudes as calculated by numerical 
method and by exact vorticity expression for an axisymmetric vortex that 
approximates the measured velocities; b = 749 mm, x^ = 7.2, 5^ = 3.0. 
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(a) Crossflow velocity. 



(b) Stream function. 

Figure 10.- Results calculated from the measured velocities; b = 749 mm 
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(c) Vorticity. 
Figure 10.- Concluded. 
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(c) Vorticity. 
Figure 11.- Concluded. 
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(a) Two wings. 
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(b) One wing, b = 749 mm. 

Figure 13.- Variation of circulation with average radius obtained from 
contours of circumferential velocity plots (not shown) . 
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(a) Crossflow velocity. 

Figure 14.- Results calculated from the measured velocities at a streamwise 
station forward of the merger point; b = 724 mm, = 1.8, = 4.2; 

"+" indicates flow visualization result for vortex center location. 
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(e) C, b = 66 cm, = 8.2 



(f) b = 66 cm, = 8.2 
Figure 15.- Concluded. 
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(f) tj;, = 3.3 
Figure 16.- Concluded 




44 


3. Recipient's Catalog No. 


1. Report No. 2. Government Accession No. 

NASA TM-78449 


4. Title and SuMitle 

laser-velocimeter surveys of merging vortices in a 

5. Report Date 

October 1978 

WIND TUNNEL - COMPLETE DATA AND ANALYSIS 

6. Performing Organization Code 

7. Author(s) 

Victor R. Corsiglia, James D. Iversen,* and 

8. Performing Organization Report No. 

A-7262 

Kenneth L. Orloff 

10. Work Unit No. 

9. Performing Organization Name and Address 

515-52-10 

Ames Research Center, NASA 
Moffett Field, Calif. 94035 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

t2. Sponsoring Agency Name and Address 

Technical Memorandum 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 

1 15. Supplementary Notes 

*Professor, Iowa State University, Ames, Iowa 

16. Abstract 


The merger of two corotating vortices was studied with a laser velocim- 
eter designed to measure the two cross-stream components of velocity. 
Measurements were made at several downstream distances in the vortex wake 
shed by two semispan wings mounted on the wind-tunnel walls. The velocity 
data provided well-defined contours of crossflow velocity, stream function, 
and vorticity for a variety of test conditions. Downstream of the merger 
point, the vorticity was found to be independent of the downstream distance 
for radii smaller than r/b = 0.05. For larger radii, the vorticity depended 
on the distance from the wing. Upstream of the merger point, a multicell 
vorticity pattern was found. 


17. Key Words (Suggested by Author(s)) 

Wind tunnel 
Vortices 

Laser velocimeter 


18. Distribution Statement 

Unlimited 

STAR Category - 02 

19. Security Qassif. (of this report) 

20. Security Classif. (of this page) 

I 21. No. of Pages I 

22. Price" 

Unclassified 

Unclassified 

I 1 

48 1 

$4.00 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 

NASA-Langley, 1978 




















National Aeronautics and 
Space Administration 

Washington, D.C. 

20546 


Official Business 

Penalty for Private Use, $300 


THIRD-CLASS BULK RATE 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 






POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 




